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Abstract 

Using the bosonization technique, we have studied a spin-1/2 magnetic im- 
purity in Heisenberg chain, and shown that the impurity specific heat and 
spin susceptibility have an anomalous temperature dependence. The tem- 
perature dependence of the impurity specific heat is: Cj m (T) ~ T M , fi = 
4 — |(1 — g) 2 , for g c < g < 1. The impurity spin susceptibility has the follow- 
ing temperature dependence: Xim(T) ~ T u , v = 3— — g) 2 , for g c < g < 1, 
and v = —1, for g < g c , where g c satisfies: 4g c = (1 — g c ) 2 , and g is a dimen- 
sionless coupling strength parameter. 
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Recently, the quantum impurity scattering of the Tomonaga-Luttinger(TL) liquid has 
been extensively studied by using different techniques [|lHT9|j- There is some controversy 



on the treatment of backward scattering of the conduction electrons on a quantum impu- 
rity or impurity-like hole in the valence band. One thinks that including backward scat- 
tering drastically changes the properties of a TL-liquid. In principle, this problem can 
be formulated by means of the Bethe-Ansatz to obtain exact results for static and ther- 
modynamic quantities. However to obtain frequency dependent quantities one has to use 
perturbative methods. Due to strong correlations between the conduction electrons and 



impurity, perturbation-theory may fail. In Ref. ||19|| , we have given a clear expression and 
some exact results within the validity of the bosonization method for a quantum impurity 
scattering of a general one-dimensional (ID) interacting electronic system, and shown that 
the system has two independent collective modes: one is drastically influenced by electron- 
electron interaction, and another is independent of the electron-electron interaction which 



is clearly observed in the numerical calculations ||18|| , while it is missed in the renormaliza- 
tion group calculations [f|]. Here we use the same method as that in Ref. [|T^] to study a 
magnetic impurity in the Heisenberg chain which can be easily treated by numerical meth- 
ods. The impurity susceptibility of this system shows an unusual temperature dependence: 
XimiT) ~ T u , v = 3 - ^(1 - g) 2 , for g c < g < 1, and v = -1, for g < g c , where g c 
is defined as that: 4g c = (1 — g c ) 2 , and g is a dimensionless coupling strength parameter. 
For an antiferromagnetic Heisenberg chain, the dimensionless coupling strength parameter 
g takes the value J2T[] g = ~, and the temperature dependence of the impurity suscepti- 



bility Xim{T) is T 5//2 . The impurity specific heat C im (T) has the temperature dependence: 
Ci m (T) ~ T M , /i = 4 — ^(1 — g) 2 , for g c < g < 1. Although the impurity specific heat as 
well as the spin susceptibility depend on the dimensionless coupling strength parameter g, 
we can still define a temperature independent Wilson-ratio in this system. 
We consider the following one-dimensional Heisenberg model 

# = [MSfsUi + sfs? +1 ) + j,s;sfa] (i) 



where Jj_ is the transverse exchange interaction strength, and J z the longitudinal exchange 
interaction strength. Using the Jordan- Wigner transformation to spinless fermion operators 
fi, = Sf+iSf = /j + exp —in J2i<i n { > S* — ft f%~\ an d after a Fourier-transformation, 
the Hamiltonian can be written in the form 

H = -J ± J2 ™s(ak)f+f k cos(ak)p(k)p(-k) (2) 

k k 

where a is the lattice constant, N is the site number, the sum over k is restricted to the first 
Brillouin zone, fk = -^J2j^ kXj fj, and p(k) is the density operator, p(k) = J2j e lkXj fj~ fj- 
One expects the asymptotic behavior of correlation functions to be determined by the low- 
lying excited states near the Fermi points at akp = ±|. Therefore one can consider a 
related model with two linear single-particle spectra tangent to the — cos(a/c) at the Fermi 
points. We may introduce fikif^k) operators to describe the fermion particles with positive 
(negative) group velocity and the associated fields, ipi(x) = ^J2k fik^ kx , i — 1,2, where 
L = aN is the length of the system. The density operator is p(k) = pi(k) + P2{k) with 
Pi{k) = J2pfik+ P fip f° r i — 1)2. For the case of spinless fermions, the J 2 -term in Eq. (0) is 
Sfc~2fc F p(k)p(— k) — > —J2k~oPi(k)P2(—k). Therefore the Hamiltonian (H) can be simplified 



as 



4 7 

H = v^Kftkfik - fM - -^£pi(*M-*) • (3) 
k iV k 

The Jordan- Wigner transformation for the spin operators on the lattice has the obvious 



generalization to the continuum situation [21 



™Y,fi + fi -™7 dy[ Pl (y) + p 2 (y) + (2a)- 1 ] = iN{ Xi ) 

l<i 

and the resulting representation for continuum spin operators is 
S-(x) = (±) 1 / 2 iMx)+Mx)]e lN{x) 

S + (x) = [S-(x)]+ (5) 
2S z (x) = pi(x) + p 2 (x) + ijjf (x)ip 2 (x) + if?2 (x)ij)i(x) 



These representations of the spin operators are different from that used in Ref . j2(J . Different 



results of Ref. JZIJ may come from the special choice for the spin operators and the Kondo 
interaction between the impurity and the conduction electrons. It is well-known that the 



bosonization representations of the fermion fields are pl|-f23 

p ±ikpx 

Mv) = -s=e^W (6) 
V 2na 

2n e -j\k\-ikx 

M*) = Ti-j- E 1 P*(*0 

k 

where the negative sign for i — 1 and the positive sign for i = 2. a^ 1 is an ultraviolet cutoff, 
which is of the order of the conduction band width. The Hamiltonian (|3|) can be written in 
terms of boson density operators 

H= 2jw L j2[ Pl (-k) Pl (k) + p 2 (-k)p 2 (k)} 

L k>0 

L k 

= - 7 2 )* J dx[{d,^ + {x)f + (d x ^(x)) 2 } (7) 

where [px(-k), Pl (k')} = [p 2 (k'), p 2 (-k)} = %5 kk r, $ ± (x) = ^(z) ± $ 2 (z), $+(x) = 
*_(*) = (g)V4$_( x)) ^$ ±(a;) = 2 7r( Pl (x)±p 2 (x)), 7 = -2J,/(*v F ). 
The boson fields $-t are conjugate variables, [(9 X $ + , = i47r5(x — y), [d x <&_(x), 

<£ + (y)] = z47r5(z-?/). 

We consider a spin-1/2 magnetic impurity scattering on the conduction fermions 

H im = j£(S + S-(0) + S-S + (0)) + j*S*S*(0) (8) 

By taking 5* + = / + , S 1- = /, S z = f + f — |, the Hamiltonian (f|) can be written as 

= 4if + (M0) + M0))e m{0) + h.c] (9) 
+Jo(f + f- i)[pi(0)+p 2 (0)+^ + (0)^ 2 (0) + ^+(0)^ 



where J5 1 = J$ /y2a and Jq = Jo/4. However, due to the appearance of the cross term 
ipi (0)i/) 2 (0) + •0 2 "(O)'0i(O) i n ©> which describes the backward scattering of the magnetic 



impurity on the conduction fermions, the problem becomes unsolvable. Therefore we have 
to find a method to eliminate this cross term. For this purpose, we define a new set of 
fermion operators 

*i(x) = -j=(Mx)+M-x)) (io) 

and the density operators Pi( 2 )(x) = \ E P Pi(2)(p)e tpx and pi( 2 )(p) = Efc ^t(2)( k + p)*i(2)(&) 
which have the standard commutation relation of the right mover, [pi(2)( — Pi(2)(9')] = 
7T$qq', [f>i(q), P 2 (g')] = 0- Now the bosonization representation of the fermion fields 



^i(2)(%) can be performed in the standard way |2l|-p3 



*^(x)= * e -*x«(*> (11) 

01(2) (X) =i T E P1(2)(P) 

where d x 4>i{2){x) = 27rpi( 2 )(x). Using Eq. (|TD|) to express the density operators pi( 2 )(a;) 

P2(a;) = 2^i(- x ) + ^(-x) - #2"(-z)# 2 (-a0 
-*^(-a;)*i(-a;)] 
Vi + (0)V 2 (0) + ^ 2 + (0)^(0) = px(0) - p 2 (0) (12) 

and defining new boson fields: 4>±(x) = <pi(x) ± 2 (x), which satisfy the commutation 
relations: [4>±{x), 4>±{y)} = iATt5{x — y), [<f> + (x), 4>-{y)} = 0, the Hamiltonians (0) and 
(U) can be rewritten as 



+ (<£'_ (x)) 2 - ^cos(0_(x))cos(0_(-x))} (13) 
cr 



ff iw = v / 2J5 L (/ + *i(0)e iiV ( ) + h.c.) 

+ 2 f(/ + /-i)(0' + (O)+0'_(O)) (14) 



where <j)±(x) = d x (f)±(x). The cross term in Eq. (|14]) has vanished, but the boson field <f>-{x) 
in (|T3|) becomes strongly self-interacting. The Hamiltonian (|14D can be rewritten as 



H im = v^J o ± (/+* 1 (0)e w (°) + /i.e.) (15) 
2tt 



+^^(/ + /->;(o) 



+ ^fe__2M (/+/ _i )0l(o) 

where for simplicity, we use J${k = 0) and J${h = 2k f) to indicate the forward and back- 
ward scattering interaction strength, respectively, because they show different effects on the 
system. We can cancel the term in ( |I5|) by the following unitary transformation 

U = exp{*(/ + / - hi^MO) + -<M0)]} (16) 

I 7T 7T 

where, 5+ = arctanf J o( fc ~°L ) g_ = arc tan(^ a %— — ), and q = (tt 1 ) 1 ^ 2 is a dimensionless 
coupling strength parameter. 

However, the backward scattering potential has a drastical influence on the fermions 
^1(2) (x), and induces the strong coupling between the fermion fields ^/i( 2 )(x) and / at the 
impurity site x — 0. In the strong coupling limit induced by the backward scattering 
potential, i.e., the phase shift £_ takes the value: 5°_ = — tt/2, and taking the gauge trans- 
formations: *i(2)(z) = ^i(2)0r)e i9l(2) , Q x - 2 = 25_(f + f - 1/2), we can rewrite the total 
Hamiltonian as 

H T = U + (H + H tm )U 

= 7T I dx{((j) + (x)) 2 + ^<j) + (x)<j) + (-x) 



in 



+ (0_(x)) H — - cos(0_(x)) cos(0_(-x))} 
or 

+ Ijj-tf+e^+tot+Wu+e^Wu + h.c.) 
Jtxol 



where <j>±{x) = fa(x) ± fa(x), d x fa( 2 )(x) = 2np m (x), pip){x) = ^i{2)( x )^m( x )- The 
last term can be easily obtained by perfoming the unitary transformation U + f + U which 
contributes the phase factor —ig5 + /(7r)<j) + (0) — i5-/(Tr)<f>-(0). In the strong coupling limit 
(51 = —tt/2), by using Eq.(ll), we can obtain the last term in (|T7|). If we redefine the 
following new fields: 
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1 



(#i(x) + * 2 (x)) 



v/2 



(18) 



"02 (x) 



1 



(#i(-z) - tf 2 (-a:)) 



where the bosonization representation of the fermion fields Vi(2)( a; ) is '01(2) C^) = 
(— L) 1 / 2 ^- 1 *!^)^ the total Hamiltonian can be written as 



where 0>(O) = (f^) 1/2 exp{-z(± + ^ ± )($ 1 (0) + $ 2 (0))}. After these transformations, the 
phase factor N(0) disappears in the total Hamiltonian (|19|). If the dimensionless coupling 
strength parameter g takes the values, g > 1, in the strong coupling limit induced by the 
forward scattering potential, the phase shift S + satisfies the relation: 5+ = —n/(2g), 0>(O) 
becomes a constant field, the total Hamiltonian ( |19"D is very similar to that in Ref. |24j] derived 
from the quantum dot. The J^-term opens a gap in the energy spectrum of the fermion /, 
therefore, in the low temperature and low energy limit, the Green's function of the fermion / 
is an exponential decaying function, the impurity susceptibility exponentially goes to zero as 
the temperature going to zero. However, if the dimensionless coupling strength parameter g 
is less than one, g < 1, the phase shift 5 + only takes the value: 5+ = —n/2, the temperature 
dependence of the impurity susceptibility has a power- law form (see below). Near the strong 
coupling critical point induced by the forward and backward scattering potentials, we have 
the leading irrelevant Hamiltonian 



H c T = V -^j ^{($' 1 (x)) 2 + 7 $' 1 (x)«l' 1 (-x) 

+ ($' 2 (x)) 2 + 7 <i;(x)$;(-x) 

+ >/2#[/V(0)+^+(0)/) 



(19) 



AH = X(f + f~- -)0' + (O) + A'(/ + /~- -)0'_(O) 



(20) 




(21) 
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Therefore, in the srong coupling limit, from fll8|), (|19D and fl2T|) , we can easily obtain the 
following correlation functions 

<0'_(O,O)0'_(O,t)>~(i)l 

< e ^ + (o,o) e -^+(o,t) > ^ (I)f (22) 



1 

For an antiferromagnetic Heisenberg chain, the exponent g equals 1/2 (7 is very large). 



<0' + (o, 0)0^0, t)>~ ( lv2 



From Eq. (|19[ ) we can obtain the Green's function of the impurity fermion 



G f it) =< f(0)f + (t) >~ < 



0, g>l 

(j) 2 -* (1 - ff)2 , 9c<9<l (23) 
e - *'*, 9 < 9c 



where e/ is the level of the impurity fermion /, and g c is defined as that: 4g c = (1 — g c ) 2 . 
The physical interpretation of the special parameter g c is that at this point the self-energy 
of the impurity fermion / induced by the interaction term in (19) has a linear frequency 
dependence. With the help of Eq . (p2|) we see that in Eq.(^) for the antiferromagnetic 
case g < 1, the A-term is dominant while for the ferromagnetic case g > 1, the A -term 
is dominant. However, in the case of g > 1, the Green's function of the fermion / is 
exponentially decaying (in Eq.([23|) we take it as zero in the long time limit), only the A-term 
be relevant. By using Eqs. (|22|), and (p3[), the correlation function of AH (Eq. (|20|) ) can be 
written as < (/+/- V 2 )(*) " (/ + /- V 2 )(0) >< t)</>' + (0, 0) >, if we omit the vacuum 

fluctuation of the fermion /, it reads 

f A 2 (V-i (H)2 , g c <g<l 
< AH(0,0)AH(0,t) >~ { t (24) 

(a 2 ^) 2 , g<g c 

while the spin susceptibility of the impurity is (omitting the vacuum fluctuation of the 
fermion /) 

< S z (0)S z (t) >~ (- t )'-^- 9)2 , g c <g<l (25) 



From Eqs. (p4j) and (0), using the relation between specific heat and free energy (Eq.(p4|) 
gives the impurity free energy), we can easily obtain the temperature dependence of the 
impurity specific heat Cj m (T) and spin susceptibility XimiT) 



which show an unusual temperature dependence in the low energy and low temperature 
limit. We can give a brief explanation about present results. It is worth noted that because 
we use an usual Kondo interaction as in [21] which is different from that in [20], we obtain 
different results from that in [20]. The symmetry of the Kondo interaction term signifi- 
cantly influences the low temperature behavior of the magnetic impurity. In generally, the 
Heisenberg chain can be described by an interacting spinless electron system. As the inter- 
actions among the electrons are repulsive, the bound state of the conduction electron and 
the impurity fermion is weakened by the repulsive interaction of the conduction electrons. 
For an enough strongly repulsive interaction g = g , where go satisfies: 3go — (1 — 9o) 2 = 0, 
the bound state at the impurity site is broken, and the impurity fermion begins to show 
a free-type behavior. For an antiferromagnetic Heisenberg chain, the dimensionless cou- 
pling strength parameter takes the value [^TJ g = |, the impurity specific heat C im (T) is 
proportional to T 7 / 2 , and the temperature dependence of the impurity spin susceptibility 
Xim(T) is Xim(T) ~ T 5 / 2 . Because in the Jordan- Wigner representation, the Heisenberg 
chain can be reduced into a very simple form, these anomalous temperature dependences of 
the impurity specific heat and spin susceptibility should be easily observed in the numerical 
calculations. Although these temperature dependences heavily rely upon the dimensionless 
coupling strength parameter g, we still have a temperature independent Wilson ratio which 
exists in the high- dimensional Kondo problems. 




(26) 




(27) 
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In summary, by using the bosonization method, we have studied in detail the low temper- 
ature physical behavior of the spin- 1/2 magnetic impurity in Heisenberg chain, and shown 
that the impurity specific heat and spin susceptibility have unusual temperature dependence 
behavior in the low energy and low temperature limit. 

The author would like to thank Prof. P. Fulde for encouragement and Dr. K. Fischer 
for valuable discussions. 
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